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ABSTRACT 



Isotropic diffusion models for Galactic cosmic ray transport put tight constraints on the maximum convection velocity in the halo. For 
a half halo height of 4 kpc the maximum convection speed is limited to 40 km/s in the halo, since otherwise the constraints from local 
secondary to primary ratios and radioactive isotopes cannot be met. The ROSAT Galactic wind observations of wind speeds up to 760 
km/s therefore constitute a problem for diffusion models. 

It is shown that such wind speeds are possible, if the diffusion coefficient in the halo is different from the diffusion coefficient in 
the disk. The radial dependence of the wind velocity was taken to be proportional to the source strength, as expected from winds 
which are sustained by cosmic ray pressure. In this case the cosmic ray density and with it the diffuse y-ray production from nuclear 
interactions are suppressed near the sources. This solves in a natural way the problem of the soft gradient in the radial dependence of 
the y-ray flux. Furthermore, the large bulge over disk ratio in positron annihilation as observed by INTEGRAL, can be explained by 
positron escape from the disk in such a model. 

Key words. Convection - cosmic rays - ISM: jets and outflows - Galaxy: halo - Galaxy: structure 



1. Introduction 

High energy cosmic rays (CRs) are thought to be accelerated 
in the shock waves of superno vae remnants (SNRs), as can 
be fo u nd in standard text books dGinzburg et al.lfl990l ; iLongaiH 
19941; ISchlickeiset] | 2002l) . but critical reviews exist (see e.g. 



Kirk & Dendvll200ll)~ Although CRs are accelerated to almost 



the speed of light they are confined to our Galaxy for about 10 7 
yrs, as is known from the presence of radi oactive nuclei w ith de- 
cay times of a few times 10 6 yrs, see e.g. dCesars kvlll980t) . Such 
"cosmic clocks" require averaged drift speeds of CRs of only 
a fraction of the speed of light. This i s possible if CRs con tinu- 
ously scatter on magnetic turbulences dGinzburg et aljl990l) and 
as a consequence perform a random walk which can be modeled 
by diffusion. Possible scattering centers are the turbulent mag- 
netic fields generated by the dilute plasma of CRs themselves. 
In the simplest case diffusion is assumed to be isotr opic with the 
same diffusion coefficient in the halo and the disk dStrong et alj 
120071) . The parameters of the diffusion equation can be obtained 
by taking a source distribution proportional to the SNR distribu- 
tion and measuring the local density and energy spectra of the 
CRs, which depend on the transport parameters, gas densities 
and magnetic fields between the source and the solar system. 

An isotropic diffusion model has been implemented in 
the publicly available GALPRO P cod43 dStrong & Moskalenkol 
Il998t iMoskalenko et alj 119981) . Details about the underlying 
propagatio n and the parameter determination can be found in 
a review bv lStrong et alj d2007l) . 

In addition to dif fusion convec tive transport modes may 
play a role for CRs dJokipiil Il976l) . Supernovae (SNs) eject 
hot gas into the interstellar medium (ISM) which can expand 
into the halo, presumably d riven by the CR pr essure from 
supernovae remnants (SNRs) dBreitschwerdtl l2QQ8h . This leads 



1 The GALPROP code 
http://galprop.stanford.edu 



publicly available under 



to a reduction of the CR density and CR interaction rate in 
the Galactic disk, most efficiently at radii with high source 
density. This mechanism has been proposed as an explanation 
for the relatively small production of diffuse y-rays near the 
sources in comparison to the y-ray production at larger radii 
dBreitschwerdt et all l2002h . known as the "soft y-ray gradient 
(SGRG) problem". 

Unfortunately, the maximum allowed convection speed in 
isotro pic propagation mo dels is restricted to a few tens of 
km/s dStrong et al.l 120071) . because otherwise the constraints 
from local secondary to primary ratios and radioactive isotopes 
cannot be met. Wind velocities as small as this will not solve 
the soft-y-ray-gradient problem. Within the isotropic diffusion 
models alternative solutions have been investigated; e.g. putting 
more sources in the outer Galaxy dStrong & Moskalenkoll 19981) 
or putting more gas there by assuming that the tracer for molec- 
ular hydrogen, the radiation from CO mol ecules, has a differen t 
proportionality factor in the outer Galaxy dStrong et al.ll2004bh . 
as motivated by the change in metallicity. A recent analysis 
of the ROSAT data on x-rays implies wind speeds of up t o 
760 km/s in the halo dEverett et al.ll2007t lBreitschwerdll2008h . 
Speeds like this are sufficient to solve the SGRG problem, but 
they would require to loosen the constraint of isotropic diffusion. 



An additional problem for isotropic diffusion models, is 
the large bulge-over-disk (B/D) ratio of the 511 MeV positron 
annihilation line, as observed by the INTEGRAL satellite 
dKnodlseder et alj 120051: IWeidenspointner et al.l 120071) . These 
low energy positrons are thought to originate predominantly 
from the decay of radioactive nuclei, produced in SNIa explo- 
sions in the bulge and in the disk. Predictions of the B/D ratio 
from the expected number of SN explosions are well b elow one, 
becau se of the high rate of SNIa explosions in the disk dPrantzosi 
120061) . However, INTEGRAL found a B/D ratio of a few. In a 
diffusion model without convection MeV positrons hardly prop- 
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agate, because diffusion is proportional to the velocity and en- 
ergy of the particle. In this case positrons annihilate close to their 
sources leading to a small B/D ratio. Convective transport in 
Galactic winds is independent of energy, so low energy positrons 
in the disk can be convected to the halo, where there are hardly 
electrons at rest to annihilate with. In this case the escape frac- 
tion of positrons from the disk can be sufficient to explain the 
observed large B/D ratio. 

In this paper an anisotropic transport model, which allows 
for ROSAT compatible convection and still meets all the con- 
straints from primary and secondary cosmic rays, cosmic clocks, 
y-rays and the INTEGRAL B/D ratio, is introduced. The model 
features different diffusion coefficients in the halo and in the 
disk. By increasing the diffusion coefficient towards the halo 
boundary a smooth transition to free escape is obtained, so the 
model becomes insensitive to the precise position of the bound- 
ary. This is in strong contrast to isotropic propagation models, 
which require a precise size of the halo for a given diffusion 
constant. 

The paper is organized as follows: in section [2] we summa- 
rize the main features of the current isotropic diffusion models 
for CR transport. Section [3] discusses the limits of the isotropic 
transport models in the light of observational constraints which 
refine our picture of CR transport. In section|4]we introduce the 
basic ideas of an anisotropic propagation model (aPM). Section 
[5] discusses the parameter determination of the aPM. The per- 
formance of the aPM with respect to charged CRs and diffuse 
Galactic y-rays is discussed in sections [6] and |7j respectively. 
Section[8]summarizes the results. 



2. Isotropic propagation models 

CRs are thought to be accelerated in shock waves of supernovae 
explosions up to velocities close to the speed of light or ener- 
gies up to 10 15 eV. However, they do not escape with such high 
speeds, but are trapped in the magnetic field of the Galaxy for a 
time of the order of 10 7 years. During this time they diffuse ran- 
domly through the Galaxy and cross the disk many times. During 
each crossing secondary CRs, like B, Be, y rays, positrons, an- 
tiprotons and ot hers can be produce d. More details can be found 
in the review bv lStrong et al.l d2007l) and references therein. 

The main features of our Galaxy are a barred central bulge 
with a diameter of a few kpc and a large spiral disk with a radius 
of about 15 kpc and an density falling exponentially in R with 
a scale length of about 2.5 kpc and in z with a scale height of 
about 0.25 kpc. Most of the gas is distributed in the disk with a 
broad maximum between R =4 kpc and 6 kpc for molecular hy- 
drogen, while the distribution of ionized hydrogen is nearly con- 
stant between R =4 kpc and 13 kpc. The supernovae remnants 
(SNR) are also distributed in the disk, but peak at a distance 
of a few kpc from the cente r with a slow fall off to larger radii 
dCase & Bhattacharvalll998l) . as will be shown later. The CRs 
form a plasma of ionized particles, in which the electric fields 
can be neglected by virtue of the high conductivity and the mag- 
netic fields form Alfven waves, i.e. traveling oscillations of ions 
and the magnetic field, where the ion mass density provides the 
inertia and the magnetic field line tension provides the restoring 
force. If the wavelength of the Alfven waves equals a multiple 
of a particle gyroradius, resonant scattering occurs, leading to 
a change in the CRs pitch angle without energy losses. Such a 
process leads to a random walk, which can be described by a 
diffusion equation (see e.g. dGinzburg et al.ll990tlLongairll994t 



ISchlickeiserj|2"002l) ) 

df p df 1 8 , df 

i = v '*w - (HV)/ + 3 V • % + ?t/ d »w (1) 

where f(r, p, t) is the CR phase space density, Dy are the com- 
ponents of the diffusion tensor for spatial diffusion, D pp the dif- 
fusion coefficient in momentum space and u is the convection 
velocity. Convective transport is possible either by a large scale 
motion of the interstellar medium with velocity V or by the ef- 
fective velocity of the Alfven waves u w . Assuming that the en- 
ergy density in Alfven waves propagating in opposite directions 
is the same u w vanishes and u — V. In addition to spatial diffu- 
sion, diffusion in momentum space can be caused by CRs scat- 
tering on moving Alfven waves, which leads to diffusive reac- 
celeration. One finds the corresponding diffusion coefficient to 
be D pp cc p 2 v 2 a ID zz , where the Alfven velocity v a is introduced 
as a characteristi c velocity of weak disturbances propagating in 
a magnetic field dStrong et alj|2007l) . 

Locally CR diffusion is highly anisotropic and occurs along 
the magnetic field lines. However, if the B-field has no preferred 
direction, i.e. if the turbulent small scale component (~ lOOpc) 
is larger than the regular large scale component, CR diffusion is 
assumed to be isotropized by pitch angle scattering on these tur- 
bulences. This is in agreement with the fact that the flux of CRs 
is the same from all directions to a high degree. 
It is usually assumed that the steady state condition is reached, 
i.e. df/dt = 0, which implies that the injection rate of CRs by 
sources equals the loss rate. CRs can be lost either by fragmenta- 
tion, decay or escape from the Galaxy. In the steady state case the 
diffusion equation for CRs can be solved numerically for given 
boundary conditions. Usually one requires the CR density to be- 
come zero above a certain halo height. 

Most primary nuclei show a power law spectrum falling with 
rigidity like £T 2 54 . This can be easily reproduced by selecting 
the injection spectrum of the primary particles, the rigidity de- 
pendence of the diffusion coefficient and the energy gains due to 
diffusive reacceleration accordingly. During their journey CRs 
may interact with e.g. the gas in the Galaxy and produce sec- 
ondary particles. This changes the ratio of secondary/primary 
particles, like the B/C ratio. From this ratio one can determine 
that the amount of matter traversed (grammage) by a CR dur- 
ing its residence time tcR- Grammage is given by p g ct CR , where 
ctcR is the path length for a particle traveling with the speed of 
light c and p g is the gas density. It was found to be of the order 
of 10 g/cirT 2 dSchlickeiseril2002t iGinzburg et al.lll990l) . which 
corresponds to a density of about 0.2 atoms/cm 3 . This is sig- 
nificantly lower than the averaged density of the disk of about 
1 atom/crrr 3 . Under the assumption of a homogeneous gas dis- 
tribution this suggests that CRs travel a large time in low density 
regions, like the halo. 

For relativistic energies the inelastic cross sections for secondary 
particle production usually do not strongly depend on the ener- 
gy of the particle. This would lead to rather flat spectra for the 
secondary/primary ratios in contrast to the observed B/C ratio, 
which shows a maximum at about 1 GeV/nucleon and decreases 
as E~ - 6 towards higher energy, as shown by the dashed line in 
the right hand side of Fig.Q] This can be accommodated by as- 
suming that energetic particles diffuse faster out of the Galaxy, 
i.e. the energy dependence of the diffusion coefficient and the en- 
ergy gains due to diffusive reacceleration are chosen accordingly 
(see the full line in the left hand side of Fig[T|). The decrease at 
low energies can be accommodated by both diffusive reaccelera- 
tion, which shifts the spectrum to higher energies, as well as con- 
vective transport mechanisms dStrong et al.l20 07). Alternatively, 
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Fig. 1. Left Demodulated local B/C ratio for different variations of an isotropic model: isotropic reference model (full line), isotropic 
model with constant (energy independent) diffusion (dashed line), isotropic model without diffusive reacceleration (dotted line). For 
energy independent diffusion (which means the escape time does not depend on energy), the B/C ratio is independent of energy 
once the nuclei become relativistic. Setting the Alfven velocity v a to zero eliminates diffusive reacceleration, which therefore shifts 
the peak in the B/C ratio to lower energies. The grey band refers to a run with Zh = 5.3 kpc. Right Demodulated local Be/ 9 Be 
ratio for different lifetimes of l0 Be: n/2 = 1.6 ■ 10 6 yrs (full line), T1/2 = 1.6 ■ 10 4 yrs (dashed line), T1/2 = 1.6 • 10 8 yrs (dotted line). 
The grey band refers to a run with Zh = 5.3 kpc. 



one could assume a strong increase of the diffusion coefficient at 
low e nergies due to damping of the Alfven waves dPtuskin et alj 
120061) . 

From the ratio of unstable/stable secondary nuclei (like 
w Be/ 9 Be) one obtains the average residence time of CRs in the 
Galaxy to be of the order of ten = 10 7 yrs dCesarsky| [T980). Fig. 
Q]shows the local 10 Be/ 9 Be-fraction assuming different lifetimes 
of w Be. The data require residence times between 1.6 • 10 6 and 
1.6 • 10 8 years. 

A widely-used program providing a numerical solution to 
the d iffusion equation is the pu blicly available GALPROP 
code dStrong & Moskalenkoll2006l) . The transport equation in 
GALPROP is not written in the form of a phase space density 
f(r, p, f), but in the form of density per unit of total particle mo- 
mentum ^(r, p, t) defined by 



V(p)/dp = 4np 2 f(p). 
The third term in Eq[T]then can be rewritten as follows 



(2) 
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Using[3]and|2]and multiplying the transport equationQ]by p 1 one 
obtains 
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(4) 

where it is assumed that the diffusion coefficient D is a scalar 
quantity with the same value everywhere and in all directions, 
i.e. the tensor D/j has only diagonal components, which are all 
equal. The full transport equation used in GALPROP can then 



be written as 
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where q(r, t) is the source term, t/ is the time scale for fragmen- 
tation and T r is the time scale for radioactive decay. 
The basic parameters for GALPROP are the injection spectrum 
parameters, the diffusion coefficient D, the convection veloc- 
ity V, Alfven velocity v a , which enters D pp , and the size of 
the transport box, needed as a boundary to solve the differen- 
tial equation. Usually one assumes the density of scattering cen- 
ters outside the diffusion box to be so small that CRs propa- 
gate freely with the speed of light to outer space (free escape). 
Consequently the CR density becomes zero at the boundary. For 
an isotropic propagation model one has to assume that the den- 
sity of scattering centers is constant inside the diffusion box and 
abruptly drops to zero outside. Of course a smooth decrease in 
the scattering centers would be a much more natural solution, 
but this requires an increase in diffusion coefficient towards the 
boundary. This is exactly what is needed in anisotropic propaga- 
tion models, as will be discussed in Sect. [4] In such anisotropic 
models the boundary of the diffusion box can be in principle at 
infinity in contrast to isotropic models, where the residence time 
increases with the size of the diffusion box for a constant diffu- 
sion coefficient. 

The propagation parameters are usually tuned to the sec- 
ondary/primary ratio and the unstable/stable ratio of locally ob- 
served charged particles. The production of secondary and ter- 
tiary particles is calculated in GALPROP using a network with 
more than 2000 cross sections. To check the validity of such a 
self-consistent propagation model at scales outside the kpc-scale 
of the collection volume of charged particles the fluxes of diffuse 
gamma rays are calculated by GALPROP using the emissivity of 
the complete diffusion box. The absorption of y-rays in the GeV 
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range is small, so information on the CR density and gas density 
even in the Galactic Center (GC) can be obtained. The present 
GALPROP models with isotropic diffusion indeed describe the 
y-ray fluxes remarkably well, if one ignores the anomalous ex- 
cess above 1 GeV, as measured with the EGRET spectrometer 
on board of the CGO dHunter et al.ll 1997b . This excess appears 
anyway to be not confirmed by the prelimi nary Fermi-LAT data 
dPorter & Collaborationl2009l:1porterll2009l) . Note, that the moti- 
vation to drop the assumption of isotropic propagation is entirely 
independent of the EGRET observations. 



3. Limits of the isotropic propagation model with 
new observational constraints 

In isotropic transport models propagation is homogeneous and 
dominated by isotropic pitch-angle scattering, in particular the 
scattering rate is assumed to be the same in the halo and in the 
disk. In the following we will discuss three of the main obser- 
vations that cannot be explained in a straightforward way by 
isotropic transport: the Galactic winds as observed by ROSAT, 
the soft y-ray gradient as observed by COS-B and EGRET and 
the large bulge/disk ratio observed by INTEGRAL. 

3.1. ROSAT: Galactic Winds 

It has been pointed out that diffusion models are incompatible 
with vertical wind velocity gradients larger than 10 km/s/kpc 
(IStrong et al.ll2.Q07T> . because for larger wind velocities the CRs 
are blown into the halo and the required ratio of times spent in 
the halo and the disc from the combined constraints from the lo- 
cal secondary to primary ratios and the radioactive cosmic clocks 
cannot be met. 

Until recently there was no direct observational evidence that 
the Milky Way's atmosphere might feature such a wind. It has 
been suggested that the Milky Way's SN rate is too small to build 
up sufficient CR pressure to overcome the Galactic gravitational 
potential. In so far, the fact that current transport models are in- 
compatible with significant convective transport modes has not 
been considered a serious deficiency. This changed in 2007 when 
it was found that the Milky Way drives a large scale wind, too: 
the ROSAT-satellite observed an enhancement of th e diffuse soft 
X-ray background emission dLevenson et al.llT997l) . This emis- 
sion can be e xplained best by a m ixed CR and thermally-driven 



wind mo del (Everett et al. 20071). T heir model is based on 



model bv lBreitschwerdt et al.l (|200"2T) . where the spatial shape of 
the wind velocity is given by the balance of the gravitational 
potential and CR pressure and basically follows the SNR dis- 
tribution. The fitted velocities range from 173 km/s in the disk 
to 760 km/s in the halo. Compared to starburst Galaxies which 
feature wind velocities up to 3000 km/s these wind velocities 
are still moderate. However, the impact of even moderate con- 
vection velocities of a few 100 km/s on CR transport is sig- 
nificant as we will discuss in subsection 15.21 Wind velocities 
which follow the SNR distribution will lead to an ^-dependent 
diffusion-convection boundary, thus significantly changing the 
R-dependence of the escape probability of CRs. In regions with 
a high density of CR sources the CR pressure is strong enough 
to overcome the gravitational forces and drives gas into the halo. 
CRs driven by this wind will leave the Galaxy earlier and conse- 
quently the interaction rate will be smaller than e.g. in the outer 
Galaxy. In this case grammage and escape time depend not only 
particle rigidity, but become a function of the Galactocentric ra- 
dius. 



Following iJokipiil d 19761) the convection-diffusion boundary 
z c can be defined by the parameter z c (R) ~ D ZZ (R, z c ) /V C (R, z c ). 
Above this boundary convection dominates and the probabil- 
ity for a CR to return to the diffusion region below z c de- 
creases exponentially. The R and z-dependence of the param- 
eters have been explicitly included. For the parameters of the 
aPM model including the ROSAT compatible convection, as de- 
termined in Sect. [5] the R-dependence of z c is plotted in Fig. [2] 
for different rigidities. Clearly, for most radii the regions of con- 
finement where diffusion dominates, i.e. below the curves, are 
much smaller than the halo boundary, set at 7.5 kpc. Convection 
was taken to be proportion al to the source distribution from 
ICase & Bhattacharvald 19981) (plotted in Fig. [2] As a result the 
region where convection dominates peaks at R k 4kpc for high 
rigidities and becomes much broader for low energies, since then 
the ratio of diffusion oc /JE 033 over convection (independent of 
energy) becomes small. Note that near the GC the convection be- 
comes small by virtue of the decrease in the source distribution 
in our case. But such a strong decrease is also expected from 
first principles, because of the strong gravitational potential in 
the GC, which will inhibit the launching of Galactic winds. 

3.2. Softy-Ray Gradients 

It is interesting to note that the model suggested by 
Breitschwerdt was not motivated by direct evidence for Galactic 
winds, but by another problem of isotropic transport models : The 
distribution of supernova remnants (SNR), which are believed 
to be the sources of CRs, peaks toward the Galactic center, as 
shown in Fig. [2] In an isotropic diffusion model the propagated 
CR distribution still strongly resembles the source distribution, 
leading to a peak in the radial distribution of diffuse y-rays, i.e. 
one observes a strong gradient. This is incompatible with the soft 
y-ray gradient as observed by COS-B and EGRElQ. A solution 
to this SGRG-problem would be to reduce the residence time 
and the corresponding CR density near the s ources, which is ex- 
actly w hat convection does, as spelled out bv lBreitschwerdt et al.l 
(|2002|) . This is indeed a more natural solution than an ad-hoc 
broadening of the sourc e distribution, shown by black full line 
in Fig. [2] as proposed bv lStrong & Moskalenkold 19981) or by in- 
creasing the gas density in the outer Galaxy by assuming a strong 
radial dependence of the prop ortionality constant X co of the CO 
tracer for molecular hydrogen dStrong et alj2004bh . Although an 
increase of Xco with Galactocentric radius is expected from the 
observed metallicity gradient, an order of magnitud e increase in 
Xco se ems large and is not supported by all data dStrong et al.l 
l2004bl) . In order to examine the impact of Galactic winds on CR 
transport, we assume a constant Xco factor as a first step and al- 
low for an increase with Galactocentric radius only if required. 



3.3. INTEGRAL: large bulge/disk ratio 

Low energy positrons can annihilate with electrons, preferen- 
tially bound to nuclei in order to prevent a large momentum 
transfer during the collision. A detailed accoun t of th e anni- 
hilation process was given by iGuessoum et alJ d2005h . Such 
positrons largely originate from the decay of radioactive nuclei 
expelled by dying stars, especially SNIa. In the case of SNIa the 
SNR core makes up a large fraction of the mass, so it has a rela- 



2 Note, that this y-ray gradient refers to a spatial feature of the diffuse 
y-rays. A possible miscalibration of the EGRET instrument will not 
modify the spatial shape of the observed radiation. 
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Fig. 2. Left Convection diffusion boundary as defined by D zz {R,z c ) = V z {R,z c ) ■ z c in an aPM for different particle rigidities. The 
full lines correspond to protons, the dotted lines correspond to electrons. For rigidities larger than 10 4 MV D zz and therefore z c is 
the same for protons and electrons. The kink at lkpc above the plane is caused by t he sudden incre ase of the diffusion coefficient in 
this model (see section|4]for details). Right Source distribution: pulsar distribution dLorimej2004l) ( blue, dot t ed), SN R distribution 
dCase & Bhattacharya 1998|) ( red, dashed) and the flattened parameterization used bv lStrong & Moskalerikoldl998l) ( black full). 



tively thin layer of ejecta, which makes it easier for the positrons 
to escape. 

Light curves, which are sustained first by the y-rays in the 
shock waves and later by the electrons and positrons, suggest 
that only a few percent of the positrons escape from the ejecta 
and can annihilate outside after thermalization. Positrons anni- 
hilating inside the ejecta will also produce y-rays, but these will 
not be visible as an annihilation line due to the successive inter- 
actions in the shock wave. 

INTEGRAL found a very strong 511 keV positron anni- 
hilation line towards the Galactic center cor responding to an 
annih ilation rate of (1.5 + 0.1) • 10 43 5 _1 dKnodlseder et al.l 
120051) . In contrast, the annihilation signal from the bulge was 
very weak. A bulge/disk (B/D) ratio of a few was quoted, 
although additional data sh ow a clear disk signal as well 
dWeidenspointner et al.l 120071) . Taking the dominant source to 
be /3 + decay of 56 Co in SNIa, one would expect a B/D ratio 
to be well below one, because of the higher mass in the disk 
and the higher rate of SNIa ex plosions expec ted in the thick 
disk as compared to the bulge dPrantzosll2006l) . An additional 
problem presents the observation of the 1.8 MeV line from 
the 26 'Al radioactive isotope, which was clearly observed in the 
bulge and th e disk by the Com ptel detector on NASA's CGRO 
observatory dDiehl et al.ll20Q6h . These nuclei are thought to be 
produced by nucleosynthesis in massive stars and yield in their 
decay on average 0.85 positrons. The observed flux of positron 
annihilation in the disk seems to be saturated already by the 
positrons from 26 A/ and 44 Ti, leaving little room for additional 
positrons from SNIa explosions in the disk. Remind that 26 Al 
has a half life of the order of 7.2 • 10 5 years, so the positrons 
from their decay are not convected away by the wind from the 
dying star. In isotropic transport models this result cannot be 
understood, because the ~ MeV positrons do not propagate 
(D zz = Drr oc j3 ■ p 6 ), so that positrons annihilate close to 
their sources in the bulge and in the disk. In order to explain 
the observed B/D ratio one has to assume an addition positron 



population. These positrons have to be confined to the bulge in 
oder to reproduce the observed B/D ratio. Recently a low-mass 
X-ray binary population showing the morphological features of 
the observed emission from t he Galactic Center has been found 
dWeidenspointner et al.ll2008h . However, the expected emission 
from these binaries is not sufficient to explain the complete 
signal from the Galactic Center and in particular it is not able 
to explain, why the emis sion observed fr om the bulge can be 
already attributed to 26 Al. iPrantzosI d2006l) pointed out that B/D 
ratios as small as 0.5 are compatible with the INTEGRAL data 
if the disk positrons diffuse sufficiently away from their sources. 
He estimated that about 50% of the ~ MeV positrons have to 
leave the confinement region below zc before slowing down. 
In section 16.11 it will be shown that transport models including 
convection can explain the large B/D ratio in a natural way. 



4. Anisotropic propagation models 

The publicly available GALPROP code was modified to allow 
for anisotropic transport. The following changes were made: 

- In GALPROP an equidistant grid is used to numerically 
solve the diffusion equation. However, for parameters vary- 
ing on small scales a fine grid is required, which would dra- 
matically increase the memory requirements and computing 
time. Therefore a non-equidistant user-defined spatial grid 
with arbitrary grid points was implemented to allow for a 
course grid spacing in the halo with simultaneously a fine 
grid in the disk. 

- The isotropic diffusion coefficient is replaced by Drr for 
transport in the R direction and D zz for transport in the z 
direction. Both diffusion coefficients may depend on spa- 
tial coordinates and may have an independent energy depen- 
dence. 

- The convection velocity depends on Galactocentric radius. 
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- To model local gradients in transport parameters and gas 
density the user can specify two regions (possibly corre- 
spond ing to the Local Bubble and the Local Fluff dFrischl 
120091) ) for which the transport parameters or gas density may 
differ from the global parameters. 

- The spatial derivatives of the propagation parameters, which 
are needed to solve the diffusion equation, have been cal- 
culated for the additional R and z dependence of the diffu- 
sion and convection parameters. The corresponding Crank- 
Nicholson coefficients have been given in the Appendices lAl 
and|B] 

The parameterization of the spatial dependence of the con- 
vection and diffusion will be discussed in the following sections. 

4.1. Parameterization for Convection 

Convection was chosen to be proportional to the /^-dependence 
of the source density to take care of the increased CR pressure 
close to the maximum of the source distribution. Note, that the 
ex act proportionality is not a p rediction of the model developed 
by iBreitschwerdt et all d2002l) . However, their predicted wind 
velocity is very similar to the source distribution. The convec- 
tion velocity is parameterized as follows: 



V c (R,z) = Q(R,0)(®(\z\-zoWo + ^-z). 

dz 



(6) 



Here Q(R, 0) is the /^-dependence of the CR source distribution, 
Zq is the initial height from where the wind is launched, © is 
the Heaviside step function, Vq is the convection velocity at z = 
Zq and dV/dz the gradient of the convection velocity. A linear 
increase in z was chosen for r easons of simplicity. The so urce 
distribution was adopted from ( ICase & Bhattacharvalll998l) : 



Q(R,z) 



R 

IO \Ro 



-p 



C-«0> — H 



(7) 



1.2, 



where Rg is the Galactocentric radius of the sun, a 
P — 3.22 and z s — 2 kpc. 

The ROSAT dat a indicate an initial velocity of 173km/s at the 
base of the wind dEverett et al.ll2007h . Therefore zo — kpc with 
a local start velocity of lOOkm/s was chosen, which results in a 
start velocity of about 170 km/s at the peak of the SNR distribu- 
tion. dV/dz was chosen to be 35 km/s/kpc leading to a maximum 
convection velocity of 591 km/s at zi, = 7.5 kpc and 733 km/s 
at 16 kpc, which is in good agreement with ROSAT. The spa- 
tial distribution of the convection velocity for these parameters 
is shown in Fig. [3] 

4.2. Parameterization for Diffusion 

With increasing distance from the plane the density of scatter- 
ing centers will decrease until resonant scattering becomes neg- 
ligible at the free escape boundary. Except for this increase in 
diffusion coefficients the spatial dependence of diffusion is not 
well constrained. Here, a linear vertical increase in diffusion 
coefficients is chosen, since the CR halo density decreases in 
first order linearly towards the halo boundary; any possible R- 
dependence is neglected. In the disk diffusion is assumed to be 
constant. In this model we use a rather large zone of constant dif- 
fusion with a half height of z ( i = 1 kpc. The diffusion coefficients 
are parameterized as 



D.„ = D 



RR 



Kfef, \z\<z d 



where po is the reference rigidity, Dq is a proportionality con- 
stant, which is treated as a free parameter, 6 is the slope of the 
power law describing the rigidity dependence of the diffusion 
coefficient and /3 = v/c is the particle velocity. The parameter 5 
is taken to be the same for all rigidities, i.e. there is no break in 
the power law of the rigidity dependence. The scattering is lo- 
cally assumed to be isotropic, i.e. D zz = Drr, but the scattering 
rate, or equivalently the mean free path, depends on the position 
in the halo. 

In this paper we assume that the regular magnetic field plays 
no role, although this is clearly a questionable assumption. In 
principle one would expect an anisotropy in diffusion coeffi- 
cients, because even in the presence of strong scattering, the 
large scale component of the Galactic magnetic field is non- 
negligible dDe Marco et al.l2007l:lc~odino & Plouinl2007l) . In the 
model presented here we choose D zz = Drr, which implies lo- 
cally an isotropic scattering, but the overall scattering rate de- 
creases towards the boundary by the positive gradient in D zz and 
D xx . For a large enough gradient in convection the convection 
will dominate over diffusion above a certain z value. The dis- 
tance from the disk above which convection dominates over dif- 
fusion can be estimated from the convection-diffusion boundary, 
which has been shown before in Fig. [2] for different rigidities 
using the parameters discussed in the following section. 



5. Parameter determination for the Anisotropic 
Propagation Model (aPM) 

In this section we discuss the parameter tuning for the aPM pro- 
posed in section [4] It is not the aim of this study to present a 
fine-tuned best fit model, which anyway would be rather short- 
lived in the light of the upcoming Fermi-LAT and PAMELA data 
releases, but rather to show that the ROSAT Galactic winds in 
principle are compatible with local CR measurements. 



Parameter | aPM | Conventional Model 


Injection Spectra 


Protons/ 
nuclei 

a\l 02/03 


1.6/1.8/2.41 
2 GeV/9 GeV 


1.98/2.42 
9 GeV 


Electrons 

P\IP\ 

In 


1.6/2.54 
4 GeV 

ection spectra ar 


1.6/2.54 
4 GeV 

; of the form 

i) 


Transport Parameters 


D 
Po 
6 

Zd 


5.3 • 10 28 
4 GeV 8 
0.33 
1 kpc 


5.8 ■ 10 28 ssi 
4 GeV * 
0.33 


V 

Zo 
dV_ 

VAR) 
v a 
Bo 


100 f 

Okpc 
35 ^f- 

s-kpc 

K Q{R) 
56*f 
6.5 /iG 


30*f 
6.1 /Ig 



Table 1. Parameters of the aPM and a conventional GALPROP 
model. 
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Convection Velocity | | Convection velocity for different height above the Galactic plane 




Fig. 3. Left Convection velocity in an aP M.V c (R,z) = Q(R,0)(V + z(d V/dz)), with V Q = 100 km/s and dV/dz = 35 km/s and 
Q(R,0) given by the SNR distribution by ICase & Bhattachary a (1998) (normalized to 1 at R ). Right Convection velocity for 
different distances from the Galactic plane. Below z=0.1 kpc only the contribution from zdV/dz plays a role. The maximum wind 
velocity was chosen to be 591 km/s (in good agreement with ROSAT). 



The optimization of the parameters follows the same path 
as for the isotropic model discussed in Sect. |4j i.e. the diffu- 
sion coefficients are chosen to best reproduce the local B/C and 
10 Be/ 9 Be ratio and the injection spectra for protons and elec- 
trons are chosen to fit the local proton and electron spectra. The 
convection velocity parameters are taken from the ROSAT data, 
as discussed above in Sect. 14.21 The most important transport 
parameters for this model can be found in Table Q] and will be 
discussed in more detail below. 



5.1. Injection spectra in an aPM 

Fig. E] shows the local and local interstellar proton (LIS) spec- 
trum, where the local one is the observed spectrum and the LIS 
spectrum is the observed one corrected for solar modulation. The 
correction fo r solar modulation was do ne in the force-field ap- 
proximation dGleeson & Axfordlll968i) . as implemented in the 
GALPLOT programB 

Fig.|4]shows the local electron spectrum. The electron injec- 
tion index was chosen to give a propagated electron spectrum 
with an index of roughly 3.3. This yield a some what softer spec- 
trum than observed by Fermi (lAbdo et a l. 2009), but agrees with 
an extrapolation of the low energy electron data. The reasoning 
behind this is that for high energies local sources may contribute 
significantly to the local electron flux, because the large elec- 
tron synchrotron losses cool electrons from distant sources effi- 
ciently. These sources are not included in our calculations, which 
means that one expects our propagated electron spectrum to be 
somewhat softer than the data indicate. 



5.2. Diffusion in an aPM 

In order to cope with the large CR residence time and the com- 
parably small amount of secondaries, CRs have to spend a cer- 
tain time in the halo. The times spent in the halo and the disk, 
denoted by f/, and f f /, respectively, are constrained by the ratio 
of secondary/primary CRs, which is most precisely measured 



for the B/C ratio and the ratio of unstable/stable CRs, which is 
most precisely measured for the 10 Be/ 9 Be ratio. Isotropic trans- 
port models can reproduce the observed B/C and 10 Be/ 9 Be ra- 
tios, because the diffusion coefficient and halo size can always 
be chosen in such a way that td / (h + td) agrees with the measure- 
ments. 

In the presence of Galactic winds the probability for a CR to 
return to the disk drops fast with increasing height above the 
disk, leading to a relative reduction of td- Thus, if diffusion is 
isotropic and homogeneous, there is a limit for the convection 
speed which still allows to reproduce the correct ratio of td and 
0,. If convection speeds larger than this value are considered a 
mechanism which efficiently returns CRs from the halo is re- 
quired. This can be a larger diffusion coefficient in the halo or a 
generally larger D zz . The gradient in D zz can be tuned to repro- 
duce the correct ratio of f/,/(f/j + td) for a given wind velocity, 
which can be seen as follows. With V c = Vq + dV/dz ■ z and Eq. 
[5] CR transport along z can be written as 



dD zz dV 

~a — V °~ ~T Z 

oz dz 



9¥ 

~z +D ~ 



d 2y 



dz \dt 



(9) 



3 The GALPLOT routine is available from 
http://www.mpe.mpg.de/~aws/propagate.html 



The first and the last term on the l.h.s. in Eq.|9]describe the drift 
of CRs due to convection and diffusion. Drift due to diffusion 
occurs when CRs are diffusing from regions with large diffusion 
coefficient to regions with small diffusion coefficient where they 
are confined more efficiently. If the sign of dD zz /ddz is chosen 
to be positive, i.e. D zz increases in the halo, then the drift due 
to convection can be compensated by the drift due to diffusion 
(at least for a certain energy range). Unlike the vertical gradi- 
ent dV/dz the initial wind velocity Vb only occurs in the first 
term. An increase in Vq can be interpreted as an decrease in 
dD zz /dz, which just corresponds to the reduction of the forward 
mean scattering length in a moving reference frame. 

Of course, in reality the situation is more complicated due to 
energy losses, the energy dependence of the diffusion coefficient 
and because dD zz /dz = for \z\ < 1, but the basic idea is that 
CRs can easily return to the disk via diffusion. 

Since diffusion is energy dependent and convection energy 
independent the return probability depends on energy. Figure 
gives an impression of the rigidity dependence of the diffusion 
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Energy/nucleon [GeV] Energy[GeV] 



Fig. 4. Left Proton flux in am aPM: local proton flux ( full line) and LIS proton fl ux ( dashed line). Right Ele ctron flux in a 
conventional aPM. Line coding as on the left. Data are taken from the CR database bv lStrong & Moskalenkol(l2009l) . 




Fig. 5. The rigidity dependence of the convection-diffusion boundary z c for R = 8.3 kpc as defined by D ;c (8.3 kpc, z e ,p) = 
V c (8.3 kpc, z c ) • Z c - The change in slope at 1 kpc corresponds to the transition from constant diffusion to z-dependent diffusion. 



convection boundary z c above which convective transport domi- 
nates. Above z = 1 kpc, where dD zz /dz > 0, the critical rigidity 
increases from 2 GeV to more than 30 GeV. This energy de- 
pendence describes quite well the B/C ratio, while keeping the 
residence time in agreement with the 10 Be/ 9 Be ratio, as shown 
in Fig. [6] Here the parameters given in Table [T] were used. The 
predicted 10 Be/ 9 Be ratio is slightly below the data indicating that 
the residence time is slightly too large, but well within the errors 
of the data. The slightly too high value of the B/C ratio below 
1 GeV CRs also indicates that low energy CRs stay slightly too 
long in the Galaxy, which again indicates that dV/dz is slightly 
too low (note that an increase in the base velocity Vo> will n °t 
improve this situation, but on the contrary, an increase in Vo will 
help to keep CRs longer in the Galaxy). A refined tuning can be 
performed as soon as the PAMELA B/C data will be released, 
since the ROSAT data allow for dV/dz up to ~ 43 km/s/kpc, 



which is somewhat higher than the value of 35 km/s/kpc from 
TableQ] 

5.3. Halo size 

As discussed before, the diffusion coefficient increases with in- 
creasing z, thus providing a natural transition to free escape inde- 
pendent of the boundary condition. E.g. increasing the boundary 
box from z=7.5 to 30 kpc does not change the residence time 
or the ratio f/,/(f/i + td)> as can be seen from the gray band on 
top of the black line for the B/C ratio in Fig. [6] which is hardly 
distinguishable from the black line for a boundary of z=7.5 kpc. 
This is in strong contrast to isotropic transport models, which are 
highly sensitive to the size of the diffusion region, as discussed 
in Sect. |2](see Fig. [TJ. Given this large difference in sensitivity 
to the boundary condition between the isotropic and anisotropic 
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B/C ratio galdef ID XX.4_N568_N594 



PHI=300MV 



Be10/Be9 ratio galdef ID XX.4_N568_N594 



PHI=300MV 
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Fig. 6 . Left Local B/C ra tio in a conventional aPM (full line) and LIS B/C ratio ( dashed line). Data: HEAO-3 dEngelmann et al.l 
Il990l) ACE (iDavisI l2000t) . Left Local beryll ium fraction in a conventional aPM. Line c oding as on the left. Data: ISOMAX 
(lHams et al.l2004l) . ACE (lYanasak et al.ll200ll) . Voyager (iLukasiakll 19991) . Ulysses (IConnellll 1998J). The narrow gray band on top of 
the black line refers to a model with an increased halo height z/,=30 kpc (see section l5~3l for a discussion ). 



models it is interesting to compare the CR density profiles for 
protons in the halo. This is done in Fig. [7] for two halo sizes in 
the aPM (zi, = 7.5 and 10 kpc) and two halo sizes in an isotropic 
model (zh = 4 and 5.3 kpc) at rigidities of 0.01 GeV and 10 GeV 
at the Sun's Galactocentric radius. Most sources are located in 
the disk, so the source density at z — kpc is highest. In the 
isotropic model low energy CRs are almost at rest and the CRs 
stay for a large fraction inside the gaseous disk (indicated by the 
vertical gray band in Fig. [7j with tails towards the halo bound- 
ary, where the density drops to zero. For the aPM this distri- 
bution is broadened at larger z-values by convection. For high 
energies the situation changes, because the diffusion starts to be- 
come more important and the diffusion-convection boundary z c 
moves to z-values of a few kpc, as indicated by the hatched area 
on the right-hand side of Fig. [7] The ratio of times spent in the 
disc and in the halo is approximately given by the area below the 
curves in the disk and the halo (a disk height of 1 kpc is cho- 
sen for illustrative purposes, note that this is not identical to the 
significantly smaller scaleheight of 250 pc). One observes that 
this changes for the isotropic model significantly, if one moves 
the boundary from 4 to 5.3 kpc, while for the aPM the changes 
are mainly in the convection zone from which only few CRs re- 
turn to the disk. Convection therefore allows us to reproduce the 
vertical distribution of the CR distribution in such a way that the 
requirements from local B/C and 10 Be/ 9 Be are met and at the 
same time reduces the model's sensitivity to the exact position of 
the boundary condition, provided that the diffusion-convection 
boundary z c is small enough. Clearly, with increasing energy the 
convection-diffusion boundary also moves further into the halo 
until it becomes comparable to the halo height Zh- With increas- 
ing energy any aPM will therefore become more and more sensi- 
tive to the position of the boundary condition. At the same time 
however, CRs become less and less confined and for very high 
energies the diffusion approximation will break down. Figure [6] 
shows that for convection speeds as expected from ROSAT a 
halo height of 30 kpc still does not yield any significant change 
on the local B/C and l0 Be/ 9 Be ratio for all energies of interest, 
i.e. the energy range in which the diffusion approximation holds. 
The radial difference between the two models is less significant, 
as shown by the two-dimensional proton distribution at 5 GeV 



in the Rz plane in Fig. [18] Both distributions are similar up to 
Galactocentric radii of 4 kpc. For larger radii the aPM falls off 
less steep than the isotropic model. The reason for this is that 
the diffusion is assumed to be independent of radius, while the 
convection decreases at larger radii, so diffusion becomes more 
important there, thus widening the distribution in R. This helps 
in solving the soft y-ray gradient problem, as will be discussed 
in more detail in section [772] 



5.4. Antiprotons and positrons in an aPM 



With the basic transport parameters fixed by the local measure- 
ments of B/C and 10 Be/ 9 Be and the CR injection spectra given 
by the local proton and electron spectra the model parameters 
are basically settled. Fine-tuning of the transport parameters can 
be done by considering additional secondary particles, like an- 
tiprotons and positrons, which originate from nucleon-nucleon 
collisions in the ISM. 

The left side of Fig. [9] shows the local antiproton and positron 
flux. Antiprotons show an excess of about 40% to 50% below 
4 GeV. This excess is ver y similar to the ex cess seen in the 
isotropic transport models dStrong et al.ll2.007b . It is possible to 
increase the local antiproton flux by increasing the local CR in- 
teraction rate via diffusion or convection. However, this would 
worsen the B/C ratio and the local positron flux. 

The right side of Fig. [9] shows the local positron flux. For 
energies below 5 GeV the model shows a slight excess in lo- 
cal positrons. This is probably the result of too efficient diffu- 
sive reacceleration driven by either a too large v a or a too small 
diffusion coefficient. On the other hand less efficient diffusive 
reacceleration would worsen the local B/C ratio. An improve- 
ment for both positrons and the local B/C ratio is expected if 
convective transport was slightly more efficient. We will see in 
subsection l7.1l that a less steep radial gradient in Galactic winds 
can help to improve local charged CRs and diffuse y-rays. 
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Proton halo distribution at 0.01 GeV Proton halo distribution at 10 GeV 




Fig. 7. The halo density profile of protons for rigidities of 0.01 GeV and 10 Ge at R„ for an aPM (black line) and an isotropic model 
(red lines). In order to allow for easy comparison, the proton density for each energy bin has been normalized to the aPM with 
Zh = 7.5 kpc in this figure. 



Proton distribution at 5 GeV 



Proton distribution at 5 GeV 





Fig. 8. The density distribution of 5 GeV protons in the Rz plane for an aPM (left) and an isotropic model (right). 
Note the different boundaries in z of 7.5 and 4 kpc, respectively. The boundary in R is 20 kpc in both models. 



6. Performance of the Anisotropic Propagation 
Model - charged CRs 

6.1. The INTEGRAL positron annihilation signal in an aPM 

As discussed in Sect. I3.3I the INTEGRAL satellite observed a 
large B/D ratio for the positron 511 keV annihilation line. In 
an aPM this is expected, because low energy particles propa- 
gate predominantly by convection and convection is large in the 
disk and small in the bulge. Positrons produced in the disk are 



transported into the halo, but stay near the sources in the bulge. 
Furthermore the bulge has a large extension in all directions, so 
even if the positrons are slowly transported by diffusion, they 
still have time to thermalize and find an electron to annihilate. In 
the disk positrons are produced predominantly in the region of 
a high source density, i.e. a region of high convection. To quan- 
tify the escape probability from the disk the SNR distribution 
given by equation|7]is taken as the source distribution for MeV 
positrons, but a somewhat larger scale height of 300pc is used, 
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Fig. 9. Left The antiproton flux in the aPM, where the solid line denote s the local flux and the dashed line the flux corre cted for 
solar modulation. Data: BESS 95-97 dOrito et al.ll2000l). CAPRI CE 98 (iBoezio et alj|200ll). M ASS91 dBasini et al.lll999t). Right 
The p ositron flux in the aPM. Data: AMS I (lAlcaraz et al.ll2000h . CAPRICE 94 (IBoezio et"alll2000t) . HEAT 94 (iDuVernois et all 
2001) 
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Table 2. Positron escape fraction for different convection velocities and diffusion coefficients. For details see text. 



which is the scale height of SNla dPran tzos 2006). The positron 
spectrum from 56 Co j8 + -decays is modeled as a rectangular func- 
tion between 0.1 MeV and 5 MeV. Table|2]shows the fraction of 
positrons above z c for different parameters of the aPM. For the 
aPM more than 89% (83.2%) of the 0.1 (4) MeV positrons es- 
cape from the Galaxy. Since CR transport in this energy range 
is mainly governed by convection, variations of 10% in the dif- 
fusion coefficient or a constant diffusion coefficient in the halo 
(aPMl-aPM3) do not change these numbers significantly. In the 
disk Vq is larger than dV/dz, so that a decrease in dV/dz by 
50% (aPM4) still yields about the same positron escape fraction. 
When decreasing Vo from 100 km/s to 30 km/s for a constant 
dV/dz = km/s/kpc (aPM5-aPM7) the positron escape frac- 
tion decreases significantly and becomes incompatible with the 
INTEGRAL requirements for V = 10 km/s (ISOl). A strong 
increase in convection of dV/dz = 100 km/s/kpc is compati- 
ble with the requirements by INTEGRAL even if the velocity at 
the base of the wind is zero (aPM8), however in this case the 
requirements by ROSAT and local charged CRs are no longer 
met. For a quasi isotropic model with convection velocity of 
dV/dz = 7 km/s/kpc and Zh — 4 kpc (ISO) the positron escape 



fraction rages from 10.7% at 0.1 MeV to only 3.1% at 4 MeV. 
Models with Vo smaller than 30 km/s are generally incompatible 
with the INTEGRAL requirements. 

In summary, models that meet the requirements by ROSAT 
automatically provide a positron escape fraction high enough 
to explain the large bulge/disk ratio observed by INTEGRAL. 
Intermediate convection velocities too small for ROSAT are al- 
ready able to meet the requirements by INTEGRAL. In the keV 
range the scale height of the propagated positron distribution is 
~ 1 — 1.5 kpc i n excellen t agree ment with the scale height of 1 
kpc adopted by IPrantzosI d2006l) in his models B and D. Note, 
that the positron escape fraction shown here is just a simple esti- 
mate. A more detailed modeling of the positron source spectrum 
would have a significant impact upon the energy dependence of 
the escape fraction. 

It should be noted that due to the smaller source strength the 
number of positrons produced in the bul ge is not suffici ent to en- 
tirely explain the signal from the bulge ( lPrantzosl2006h . Even in 
an aPM an additional process or an additional source population 
in the bulge is required to explain the observed emission from the 
bulge. IPrantzosI d2006l) suggested that a fraction of the positrons 
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escaping the disk may even be channeled by the poloidal field 
to the bulge where the positrons then would be confined by the 
large magnetic field. A detailed study of the efficiency of such a 
channeling process would require to adapt the anisotropy in the 
diffusion coefficient according to the direction of the magnetic 
field in the halo and take care of the positron confinement in the 
bulge by a decreased diffusion coefficient. In this rather quali- 
tative approach we refrain from any fine-tuning of the diffusion 
coefficients. 

Recently it has been found that an additional low mass x-ray 
binary (LMXB) population seems to reside in the bulge region 
showing even the morphological features of t he observed an- 
nihilation signal dWeidenspointner et al.l 12008b . The additional 
positrons expected from the LMXB are not sufficient to account 
for the complete emission observed from the bulge and this also 
cannot explain why there is almost no annihilation signal from 
positrons from 56 Co from the disk, but added to an aPM this ad- 
ditional LMXB population might nicely explain the B/D ratio. 
With wind velocities taken to be simply proportional to the 
source distribution we already find e xcellent agreement with the 
model presented by iPrantzosI d2006l) . The problem of the large 
B/D ratio for positron annihilation is thus intimately related to 
the propagation of positrons. 

7. Performance of the Anisotropic Propagation 
Model - y-rays and radio emission 

7.1. Diffuse y-rays 

In the light of the upcoming release of the Fermi-LAT data on 
diffuse y-rays a cross-check with diffuse y-rays is difficult. The 
EGRET data do not agree with t he preliminary Fermi-LAT dat a 
even in the lowest energy range dPorter & Collaborationll2009b . 
In the absence of other options we continue to use the EGRET 
data as a cross-check, keeping in mind that the softer Fermi 
data will probably require a larger contribution from inverse 
Compton (IC) and bremsstrahlung. 

The diffuse Galactic y-rays are shown in Fig. fTOl the regions 
used in this figure are the regions as introduced in dstrong et alj 
I2005I) . The latitude and longitude ranges are given in Table [3] 
The y-ray prediction of the aPM in general is similar to the pre- 
diction of an isotropic model with the same source distribution 
and constant Xco (shown in Fig.[10]as the dotted blue line). For 
region A (GC) the flux in the isotropic model is somewhat larger 
than for the aPM. This is the result of a too high proton density 
close to the sources in the absence of convection. For the other 
regions the y-ray flux in the aPM is insignificantly larger than in 
the isotropic model, which is mainly the result of slight devia- 
tions in the propagated proton and electron spectra. The latitude 
and longitude profiles for the inner Galaxy are presented in Figs. 
E]and|T2]for 100 MeV < E < 500 MeV. The model shows the 
same deficiency of IC emission at intermediate latitudes as as the 



Region 


Longitude Z 


Latitude \b\ 


Description 


A 


330-30 


0-5 


Inner Galaxy 


B 


30-330 


0-5 


Disk without inner Galaxy 


C 


90-270 


0-10 


Outer Galaxy 


D 


0-360 


10-20 


Low longitude 


E 


0-360 


20-60 


High longitude 


F 


0-360 


60-90 


Galactic Poles 



Table 3. The longitude and latitude of the six sky regions shown 
in Fig. [TO] 



conventional isotropic model, this can also be seen from region 
D in Fig. [10] Note, that we use a simple power-law extragalactic 
background model. A more detailed model may yield a better 
description of the y-ray flux from the halo and at intermediate 
latitudes. However, this extragalactic background model should 
be determined from the Fermi-LAT data. The longitude distri- 
bution at intermediate latitudes, shown in Fig. Qj] reveals that 
the deficiency almost exclusively results from a lack of diffuse 
y-rays from the GC region for both, the aPM and the isotropic 
model. Since bremsstrahlung and emission from 7r°-decay are 
rather flat for intermediate latitudes, an increased contribution 
from IC would improve the model prediction. Independent of 
direction the model prediction below 100 MeV is somewhat too 
low, which also indicates that the contribution from IC is under- 
estimated. The deficiency in IC can be remedied by assuming 
that the Galactic electron density is somewhat larger than the lo- 
cal electron density, thus leading to a larger contr ibution from 
IC at lo w energies, as has been done previously bv lStrong et al.l 
d2004al) . Figure [14] shows the y-ray spectra for a Galactic elec- 
tron density increased by a factor 1.5, while the Galactic proton 
density remains unchanged. Figures [15] and [T_6] show the longi- 
tude and latitude profiles for the inner Galaxy for an aPM with 
the Galactic electron density increased by a factor 1.5 and the 
proton density kept constant. The right side of Fig.[l3]shows the 
longitude profile at between 150 and 300 MeV for intermediate 
latitudes. 

Whether or not the Galactic electron density can be different 
from the local electron density strongly depends on the trans- 
port parameters and the electron energy losses. Another possi- 
bility would be an untraced compo nent of non-equilibrium gas 
dBreitschwerdt & de Avillea 120061) . In this case an additional 
bremsstrahlung contribution from electrons in the hot gas would 
be expected. This contribution would be most significant at large 
latitudes in the source region, where the relative component of 
the hot gas blown out by SNs is expected to dominate. This is 
exactly what is required by the latitude and longitude profiles in 
Figs.[ITJand[T3] 

Alternatively an increase in the proton density above the plane 
in the inner Galaxy would lead to a more pronounced peak in the 
7T° emission at intermediate latitudes. This could be achieved by 
a smaller diffusion coefficient in this region which would be well 
motivated by the decrease of the regular magnetic field above the 
disk. However, this will not improve the model prediction below 
100 MeV, so that additionally a softer electron spectrum would 
be required. 

As discussed earlier the aPM is almost completely independent 
of the position of the boudary condition, i.e. the halo size, pro- 
vided, that the boundary is positioned well outside the diffusion 
zone limited by z c - In this model it is therefore possible to ex- 
amine the impact of additional IC emission from an extended 
halo with z;,=100 kpc without loosening the constraints on local 
CRs. The gain in additional IC at intermediate latitudes is at the 
% level and thus negligible in the total flux. An increased halo 
therefore will not improve the y-ray emission from the halo sig- 
nificantly. An additional untraced gas component in the source 
region, a larger Galactic electron density or possibly an addi- 
tional electron population appears to be required both in the 
isotropic model and the aPM. 

7.2. Soft y-ray gradient 

The aPM, as well as the isotropic models, require a somewhat 
increased Galactic electron density. In the isotropic case this 
leads to a pronounced peak in the Galactic IC emission from 
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Fig. 10. Diffuse y-rays for the six different sky regions as defined in IStrong et al.l d2004ah . Line coding: bremsstrahlung (light 
blue dashed), inverse Compton (green long dashed), 7r -decay (red l ong dashed-dotted), total (blue full). The pink full line is the 
extragalactic background model according to lSreekumar et all (Il998l) . The EGRET data are corrected for the point-spread function 
(PSF), for regions A, B and C the uncorrected EGRET data are shown as the grey band. Also shown is the total y-ray flux as 
predicted by an isotropic model with Do = 5.8 • 10 28 cm 2 /s at 4 GeV and v a = 30km/s (dotted blue line). 
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Galactic latitude 
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Fig. 11. Latitude profiles for the inner Galaxy (|/| < 30.5) for the aPM compared to the EGRET data between 100 and 500 MeV.The 
EGRET-excess above 500 MeV is not confirmed by preliminary FERMI data, so only data below 500 MeV are considered here. 
Line coding as in FigfTOl 




Fig. 12. Longitude profiles for the Galactic disk (\b\ < 5.5) for the aPM compared to the EGRET data between 100 and 500 MeV. 
Line coding as in FigfTOl 



X1Q" 6 galdeMD XX A N568 x10" 6 galdel ID XX.4 N644 x1 0" 6 galdel ID XX.4_N56B (electrons scaled) 




Galactic longitude Galactic longitude Galactic longitude 



Fig. 13. Longitude profiles for region D in the aPM (left), an isotropic model (middle) and the aPM with the Galactic electron 
density increased by a factor 1 .5 (right) for y-rays between 150 MeV and 300 MeV. The full blue line is the sum of the contributions 
from inverse Compton (green dashed), bremsstrahlung (light blue dotted), 7r°-decay (red fine-dotted) and extragalactic background 
(purple full). 



the GC, as can be seen from the top row in Fig.[l7] As a possible 
solution a significantly flatter source d istribution has been pro- 
posed by IStrong & Moskalenkol (Il998l) . shown in Fig. [2] as the 
black full line, which reduces the contribution from 7r°-decay 
and bremsstrahlung in the GC. However, this source distribu- 
tion has to be chosen ad hoc according to what is expected from 
y-rays. An increase in the Xco scaling factor towards the outer 
Galaxy is a more likely explanation, but the gradients required 
by the EGRET data are on the limit of what is expected from the 
increase in metallicity dStrong et al.ll2004bl) . 



As mentione d previously the self - consis tent Galactic wind cal- 
culations by iBreitschwerdt et all d2002l) predict a softer y-ray 
gradient, because the CR escape time varies with Galactocentric 
radius depending on the local source strength. Enhanced particle 
injection by the sources therefore results in enhanced CR escape 
and thus smoothens the propagated CR distribution, y-rays from 
7ro-decay predominantly originate from GeV protons. Figure [T8l 
shows the radial distribution for 1-5 GeV protons in an aPM and 
an isotropic model with Zh - 4 kpc, both models use the SNR 
distribution as the source distribution. In the isotropic model the 
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Fig. 14. Diffuse y-rays for the six different sky regions as defined in IStrong et alJ d2004al) for an aPM with the Galactic electron 
density increased by a factor 1 .5. Line coding: bremsstrahlung ( light blue dashed), inverse Compton ( green long d ashed), 7r°-decay 
( red l ong dashed-dotted), total ( blue full).The pink full line is the extragalactic background model according to ISreekumar et all 
dl998l) . The EGRET data are corrected for the PSF. 
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Fig. 15. Latitude profiles for the inner Galaxy (|Z| < 30.5) for the aPM with the Galactic electron density increased by a factor 1.5 
compared to EGRET data between 100 and 500 MeV. Line coding as in FigfPfl 
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Fig. 16. Longitude profiles for the Galactic disk (\b\ < 5.5) for the aPM with the Galactic electron density increased by a factor 1.5 
compared to EGRET data between 100 and 500 MeV. Line coding as in FigfPfl 



propagated proton distribution still resembles the strong peak of 
the source distribution, which will lead to problems in the y- 
ray production rate unless a strong increase in Xqo is assumed. 
Protons in an aPM are significantly flatter than protons in the 
isotropic model and a constant Xqo scaling factor is almost con- 
sistent with the observed flat y-ray gradient, even if the Galactic 
electron density is increased to match the IC emission from inter- 
mediate latitudes. The bottomrow of Fig.fTTlshows the longitude 
profile for the Galactic disk in this case. The improvement due to 
the flatter profile of emission from 7r°-decay and bremsstrahlung 
is clearly visible, but above 100 MeV a slight excess in emission 
from 7r°-decay and bremsstrahlung is still visible. A fine-tuned 
radial dependence of the convection velocity or a rather soft gra- 
dient in Xco appears to be compatible with the data. 

7.3. Radio emission in an aPM 

We checked that the electron energy losses via synchrotron ra- 
diation are reasonably well described by calculating the radio 
emission in longit ude an latitude and comparison to the 408 
MHz data fro m thelMaslam et al.l (1 19821) sky map. 

Following lMoskalenko et all ( 1 19981) we use the following pa- 
rameterization of the total regular magnetic field for the calcula- 
tion of the electron energy losses: 

B reg (R) = B exp C R ~ Ro ]exp[-^], (10) 

«B ZB 

with Rb = 10 kpc and zb - 0.2 kpc. We choose Bp - 6.5 /iG 
in order to best reproduce the lHaslam et al.l (1 19821) all sky map. 
Figure[19]shows the latitude and longitude profile of synchrotron 



radiation in an aPM at 408 MHz and the synchrotron spectrum. 



8. Conclusion 

A new model for Galactic cosmic ray transport is presented, 
which allows for significant convective transport as expected 
from the Galactic winds deduced from the x-ray data from the 
ROSAT satellite. The model has been realized by modifying the 
publicly available GALPROP code, which up to now allowed 
only isotropic transport. Isotropic transport models, which fea- 
ture globally constant transport parameters, can only accommo- 
date negligible convection speeds, since with convection the CRs 
are driven away from the disc, thus not returning often enough 
to the disk to produce secondary particles. 

A possible way to allow for large convection is to allow the 
diffusion in the halo to be different from the diffusion in the disk, 
i.e. to give up isotropic propagation with globally constant trans- 
port parameters. The GALPROP code was modified in the fol- 
lowing way: 

- the Galactic winds were assumed to be proportional to the 
CR source distribution, which was taken to be the SNR dis- 
tribution 

- the mean free path of CRs - and therefore the diffusion coef- 
ficient - in the halo was assumed to increase linear with the 
distance from the disk. Such a choice is motivated if the scat- 
tering centers are produced by the CRs themselves (via the 
irregular magnetic field components in the plasma) and the 
density of CRs drops approximately linearly in the halo. 
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Fig. 17. Longitude profile for the disk (1.5 > \b\) for an isotropic model (top row) and an aPM (bottom row). For both models the 
Galactic electron density is increased by a factor 1.5 to meet the requirements from y-rays at intermediate latitudes. 



30 




R[kpc] 

Fig. 18. Radial distribution of GeV protons for the aPM and an isotropic model. The source distribution is in both cases that of 
ICase & Bhattacharvaldl998l) . 



Fixing the magnitude of the convection speed to the wind speeds 
suggested by the ROSAT data, the increase in diffusion coeffi- 
cient in the halo can be fitted from the amount of secondary pro- 
duction (from B/C ratio) and the residence time of CRs (from the 
cosmic clocks, in this case the 10 Be/ 9 Be ratio). It is shown that 
such a model is consistent with all available CR data, including 
not only the ROSAT data on convective winds, but also the large 
bulge/disk ratio of the positron annihilation line as observed by 
the INTEGRAL satellite. In isotropic propagation models one 
has to introduce a new source of MeV positrons in the bulge, 
like e.g. dark matter annihilation or tunneling the positrons from 
the disk to the bulge via regular magnetic fields in the halo. In 
the anisotropic model presented here the large B/D ratio is nat- 
urally explained by the energy independent convective transport 
of low energy positrons from the disk to the halo, where there are 



no electrons to annihilate with. Convective transport is absent in 
the bulge, because the gravitational potential is too strong there 
to launch Galactic winds. 

An additional interesting feature of the present model is 
the smooth transition to free escape of CRs, because of the in- 
crease in mean free path with increasing distances from the disk. 
Therefore the boundary condition can be moved to infinity in 
contrast to isotropic propagation models, where the boundary 
condition is fine-tuned to get the correct residence time of CRs 
inside the Galaxy. 
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Fig. 19. Synchrotron latitude profile (top) and synchrotron lon- 
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for z < 0. 

For transport in momentum space the coefficients read 



111 

At 



^ =0 



02 = V(Rj,z i+ i)-V(Rj,Zi) pi 

Af z i+ \ - Zi 3(p, + i - pj) 

a 3 V(Rj,z M )-V(Rj,Zi) p M 



(A.3) 



At 



Zi+l - Zi 



3(/?;+i - Pi) 



for z > and 
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02 = V(R h Zi-i)-V(Rj,Zj) P i 

At Zi - zi-i 3(pi+i - pd 

a 3 V(Rj,Zi-i)-V(Rj,Zi) p M 



(A.4) 
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for z < 0. 
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Appendix B: Crank-Nicholson coefficients for 
anisotropic diffusion 

Follo wing the notation of the GA LPROP explanatory supple- 
ment dStrong & Moskalenko 2006) we find the Crank-Nicholson 
coefficients for the R and z dependent diffusion coefficients 
D RR {R, z, p) and D U (R, z, p) to be 

ai_ _ D RR (Ri,Zj,Pj) D RR (Rj,Zj,Pj) 
At ~ (R M - Ri)(Ri - R~0 ~ Ri(Ri-Ri-i) ~ 
_ D RR (R i+ u zj,pj)- D RR (Ri- \,Zj,Pj) 

(R i+ i -Ri-i) 2 ' 

02 _ D RR (Rj,Zj,Pj) D RR (Rj, zj, Pj) 

At ~ (R i+l - Ri)(R i+l - Rd + (R i+l -Ri)(Ri-Ri^y 



(B.l) 



(LAMBDA). Support for LAMBDA is provided by the NASA Office of Space 
Science. 



Appendix A: Crank-Nicholson coefficients for 
R-dependent convection 

The propagation method use d in GALPROP can be found in 
Strong & Moskalenko ( 1998). For the numerical solution of the 



transport equation the Crank-Nicholson implicit method is used At (z,i+\ — zi)(Zi - Zi-i) 
dPress etal.lll992l) . Fol lowing the notation of the GA LPROP ex 



»3 _ D RR (Rj+i,Zj,Pj) - D RR (Ri-uZj,Pj) _ 
AT ~ (R i+l -Ri-0 2 

D RR (Rj,Zj,Pj) + D RR {R i+ i,Zj,Pj) 

Ri(Ri+i - Ri-i) (Ri+i - Ri) 2 
for transport in R direction and 

a 1 _ DjRi. zi, pj) D zz (Rj, z i+ u Pj) - l>, {R,. Zf _i) 



(Zi+l - Zi-i) 2 



planatory supplement dStrong & Moskalenko 20Q6|) we find the 
Crank-Nicholson coefficients for R-dependent convection in z- 
direction to be 



(B.2) 



Ol = V(Rj,Zj-l) 02 = V(Rj,Zj) as = 

Af zi - Zi-i ' Af Zi -Zi-\ ' At 



for z > and 



(li an 
— =0, — 
Af At 



V(Rj,zd _a 3 
Zi+i - Zi ' At 



V(Rj,Zi+i) 

Zi+l - Zi 



(A.l) 



(A.2) 



02 = D zz (Rj,Zj,Pj) _ D zz (Rj,Zi, P j) 
At (zi+i - zd 2 (zi+i - Zi)(zi - zi-i ) ' 

03 _ Dzz(Rj,Zi+i,Pj) - D zz (Rj,Zi-i,Pj) + D :z (Rj,Zi,pj) 
Af (zi+i-Zi-i) 2 (zi+i-Zi) 2 

for transport in z direction. In the limit of an equidistant grid and 
constant diffu sion our Crank-Nicholson co efficients agree with 
those used in dStrong & Moskalen ko 2006) except for a factor 2 
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in ct\ and for transport along R. However, deriving these coef- 
ficients for the case of an equidistant grid and constant diffusion 
from their Eq. 25 we find 

a\ _ Ri-AR 

At ~ DrR WaR?' 

and 

«3 Rj + AR 

At RR Ri(AR) 2 

in agreement with our coefficients. The coefficients for diffusive 
reacceleration remain unchanged, because here only derivatives 
in momentum space occur. However, if D z: > Drr one has to 
keep in mind that v a has to be considered an effective parameter, 
scaled with the anisotropy in diffusion, e.g. v a — v a ■ -^dZ/Drr 
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